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Abstract 

We consider the motion of a compressible, viscous, and heat conducting fluid in the regime of 
small viscosity and heat conductivity. It is shown that weak solutions of the associated Navier- 
Stokes-Fourier system converge to a (strong) solution of the Euler system on its life span. The 
problem is studied in a bounded domain R C on the boundary of which the velocity held 
satisfies the complete slip boundary conditions. 
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1 Introduction 

A formal limit of vanishing viscosity and heat conductivity in the fluid models based on the principles 
of continuum mechanics gives rise to the Euler system 



dtQ + div^{gu) = 0, 

(1.1) 


+ div,j,(pu (g) u) + VxPm{q, '&) = 0, 

(1.2) 

dt + QeM{Q,'&) 

J+div^, Qp|up + peM(p,'d)^ u + pm(2,'(9)u =0, 

(1.3) 
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describing the time evolution of the basic macroscopic quantities: the mass density g = g{t,x), the 
absolute temperature -d = 'd(t, x), and the velocity held u = u(t,x). The symbols pm, denote 
the (molecular) pressure and the associated specihc internal energy, respectively. Solutions of system 
fll.llll. 3 p are known to develop singularities in a hnite time lap even if the initial state 

£*(0, •) = Po, ^9(0, •) = 'do, u(0, •) = uo (1.4) 

is regular, represented by smooth functions. Accordingly, solutions of fll.llll.2p are usually under¬ 
stood in the weak sense, where all derivatives are interpreted as mathematical distributions. 

Unfortunately, the class of weak solutions is too large to secure uniqueness and/or continuous 
dependence of solutions on the data. To remedy this problem, several additional ad'missibility criteria 
have been proposed, among which the entropy inequality 


+ diV3.(^SM(^, »^)u) > 0, 


( 1 . 5 ) 


where sm 


SM^g, d) is the specihc entropy related to cm and pm through Gibbs’ equation 


dDsM{g,d) 


DeM{g,d) +PM{g, d)D 



( 1 . 6 ) 


Although the entropy and similar admissibility criteria based on the Second law of thermodynamics 
have been partially successful when applied to problems in the 1-D geometry, see Bressan [2], Dafer- 
mos [5], they failed in identifying the relevant solution in the natural 3-D setting, see DeLellis and 
Szekelyhidi [7], Chiodaroli, DeLellis, and Kreml [3], Chiodaroli and Kreml [1] 

In the light of the above arguments and in accordance with the general approach advocated by 
Bardos et ah [1], the physically relevant solutions to the Euler system fll.lH1.3p should be identihed as 
asymptotic limits of solutions to more complex problems - primitive systems - describing the evolution 
of “real” huids. As an example of such primitive problem, we consider the Navier-Stokes-Fourier 
system in the form 

divj,(pu) = 0, (1.7) 

dt{gu) -F div^(pu <8 u) -h VxP{g, d) = div^§(t9, V^^u) - Au (1.8) 

dt{gs{g, d)) + div^(ps(p, d)u) + div^ = a, ^ V^^u) - ^ , (1.9) 

where S{d, V 3 ;U) is the viscous stress tensor given by Newton’s law 


S(?9, V^u) = i' 


fi{d) ( Va;U -F V^u - -div^ul) -h p{d)div,^ul 


, p>0, 


( 1 . 10 ) 
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and q = q('d, Vx'>^) is the heat flux determined by Fourier’s law 

q = uj > 0. 


( 1 . 11 ) 


The pressure p = pm + Pr is augmented by a radiation component, 

p{g, d) = pm(p, '&) + Pr{q, "&), Pr{q, '&) = a > 0, (1.12) 

while the internal energy reads 

^4 

e{g,i!}) = eMig,'!^) + eR{g,'d), cr = a —. (1.13) 

Q 

Accordingly, the speciflc entropy is 

An 

s(p, '&) = SM{g, 1^) + SR{g, ^9), SR{g, ^9) = SR{g, 79) = —(1.14) 

Sg 

see P, Chapter 1] for the physical background. 

Our main goal will be to show that solutions of the Navier-Stokes-Fourier system fll.7H1.9l) con¬ 
verge to those of the Euler system fll.lH1.3l) in the asymptotic regime 

z/, cn, a, A —)■ 0, (1-15) 

on condition that the limit solution is smooth. In order to avoid the well-known and still unsur- 
mountable difficulties connected with the presence of a boundary layer (see e.g. the survey by E [8]), 
we impose the complete slip boundary conditions for the velocity 


u • n\on = 0, [S(79, V^u) ■ n] x n|an = 0, 


( 1 . 16 ) 


accompanied with the no-flux condition 

q(79, Vx^9) ■ n|aQ = 0, (1.17) 

where n denotes the outer normal vector to dQ. 

Besides a vast amount of literature concerning the inviscid limit for the incompressible Navier- 
Stokes system, see e.g. Kato [15], Temam, Wang [21], [22], Wang, Xin, and Zang [23], the survey 
articles by E [8] , Masmoudi [18] , and the references cited therein, much less seems to be known in the 
context of compressible fluids. There are results on stability of certain special solutions like shock or 
rarefaction waves, see Hoff and Liu ]12], Hoff and Zumbrun ]T3], Li and Wang [16] as well as studies 
of the related linearized problems, Xin and Yanagisawa [2S]- Probably the closest result to ours has 
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been recently obtained by Sueur [20] for the barotropic Navier-Stokes system, related issues were 
discussed by Wang and Williams ^M\ . 

Our approach is based on the relative energy (see Dafermos [6] for a closely related concept of 
relative entropy) and the associated concept of dissipative solution for the full Navier-Stokes-Fourier 
system introduced in HD]. A similar strategy has been already used in m in the context of inviscid- 
incompressible limits, and in HDl, were the property of weak-strong uniqueness was established. In 
comparison with these problems, the purely inviscid, or, more precisely, zero dissipation limit features 
new additional difficulties: 

• In the inviscid-incompressible limit studied in nn, the initial data are prepared, in particular, 
the density and the temperature in the primitive system are a priori known to be close to their 
limit (constant) values. Such a piece of information is not available for the vanishing dissipation 
limit. 

• In contrast with the situation in [10], where weak and strong solutions of the same (viscous) 
problem are being compared, the uniform bounds based on the presence of viscosity and heat 
conductivity are lost in the vanishing dissipation limit. 

As a consequence of the afore mentioned difficulties, our result is path dependent - the rates of 
convergence to zero of the singular parameters n, uj, a, and A are interrelated in a special way specified 
in Section |3l For the same reason, the presence of the “damping” term — Au in the primitive system 
is necessary to control the amplitude of certain quantities, notably the velocity, on the hypothetical 
vacuum zones created by vanishing density. 

The paper is organized as follows. In Section [21 we summarize the necessary preliminary material 
including the main hypotheses and the concept of weak solution for the Navier-Stokes-Fourier system. 
In Section [3l we state our main result. The relative functional together with the associated relative 
energy inequality adapted to the present setting are introduced in Section 01 The available uniform 
bounds on the family of solutions to the primitive system are collected in Section O The proof of 
convergence is completed in Section [6l Possible extensions are discussed in Section [71 

2 Preliminaries, weak solutions 

We adopt the concept of weak solution to the Navier-Stokes-Fourier system fll.7Hl.ll]) . fll.lbp . fll.l7p 
introduced in [9l Chapter 2]. More specihcally, 

• the equations fll.7p . fll.Sp . together with the boundary conditions fll.lbp . are understood in the 
sense of distributions; 
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the entropy balance fll.91) . with fll.lTp . is relaxed to an inequality 


dt{Qs{Q, t9)) + div^(ps(p, 79 )u) + div^ ( = a, ^ V^u) - ^ 


satished in the sense of distribution; 
the system is augmented by the total energy balance 
■1 


p|u| + qe(g, d) 


Jq 12 
for a.a. r G [0, T], 


(r, ■) dx + X [ [ |up dx < [ 
Jo Jo. Jn 


:^o|uoP + Poe(po) "*^0) 


( 2 . 1 ) 


dx (2.2) 


Remark 2.1 As a matter of fact, the original definition in ja Chapter 2] does not include the 
damping term —Au, and, accordingly, stipulates equality rather than inequality in the total energy 
balance In view of the anticipated lack of compactness of the velocity field on the vacuum, 

however, the inequality seems more adequate in the present setting. 


The existence theory developed in [HI Chapter 3] requires certain structural restrictions to be 
imposed on the constitutive relations listed below: 


The pressure p takes the form fll.121) . with 


p^(^,^)=^5/2p 


J ’ 

where P G (^^[0, 00 ) fl (7^(0, 00 ) satishes 

P(0) = 0, P'{Z) > 0 for all Z > 0, 


0 < 


inz)-p'(z)z 


= Ron > 0. 


In agreement with Gibbs’ relation fll.6l) . we take 




and 


2 ” p J 

Q 


Sm{s,'^) = S 


,^ 3/2 


(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 
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where 


( 2 . 8 ) 


3 jP(Z)-P'(Z)Z 
' ^ 2 

In addition, we require 

lim S{Z) = 0. (2.9) 

Z —^CXD 

• The viscosity coefficients in (ll.lOh are continuously differentiable functions of e [0, oo) satis¬ 
fying 

|/i'(d)| < c, //(I -I- "d) < /u('d), 0 < ? 7 (d) < 77(1 -|- d) for all d > 0 (2-10) 

for certain constants /u > 0, 77 > 0. 

• The heat conductivity coefficient in Fourier’s law (ll.llh satisfies 

K e cxo), k( 1 -|- d^) < K(d) < F(1 -|- d^) for all d > 0 (2-11) 

for certain constants ft > 0, k > 0. 

The interested reader may consult 0 Chapter 2] for the physical background as well as possible 
relaxation of these restrictions. The important fact is that imposing the hypotheses fl2.3112.9p . to¬ 
gether with fll.l21irT4ll . render the Navier-Stokes-Fourier system mathematically tractable, at least 
in the context of weak solutions. Specifically, as shown in 0 Chapter 3, Theorem 3.1], the problem 
(II. 7111.171) admits a global-in-time weak solution for any choice of finite energy initial data [po, do, uq] 
satisfying the natural constraint po > 0, do > 0. 

Remark 2.2 It is worth noting that the problem considered in Chapter 3] does not involve the 
damping term — Au in the momentum eguation. However, the existence proof can he easily adapted 
to the present setting yielding the total energy ineguality in place of the eguality asserted by J2l 
Theorem 3.1]. 


3 Main result 


Following [To], we introduce the ballistic free energy 

Heig, d) = p(e(p, d) - 0s(p, d)), 
and the relative energy functional 


S (p,d, u|r, 0,U) = f 

J 


-plu 

2^' 


U|2 + 7/e(2,d) 


dHe{r,Q) 

dg 


(p-r) -i7e(r,0) 


dx. 


(3,1) 
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As shown in 0 Chapter 5, Lemma 5.1], the functional £ can be viewed as a kind of “distance” 
between the quantities [p, u] and [r, 0,U]. Specihcally, for any compact set K C (0, oo)^, there 
exists a positive constant c{K), depending solely on the structural properties of the thermodynamic 
functions stated in fl2.3112.9l) such that 


£ 'd, u r, 0, > c{K) 


|p — rp + I'd — 0p + |u — Up if [p, •d] e K, [r, 0] G iC 


1 + p|u - Up + pe(p, ■d) + p|s(p, •d)| if [p, •d] G (0, co)^ \ K, [r, 0] G K. 

(3.2) 

In order to state our main result, we need strong solutions to the Euler system (ll.im.3l) supple¬ 
mented with the boundary condition 

u-11190 = 0. (3.3) 

Note that the total energy balance (ll.3|l can be equivalently reformulated as the entropy balance 
equation 

dtigsM^Q, ^)) + div^igsMig, ^)u) = 0, (3.4) 

or the thermal energy balance 


Cv{g, ■&) {dtigd) + dwo^ig'&vL)) + div,^u = 0, with c^(p, , (3.5) 

as long as the solution of the Euler system remains smooth. 

A suitable existence result for the Euler system with the slip boundary condition fl3.3p was 
obtained by Schochet m Theorem 1]. It asserts the local-in-time existence of a classical solution 
[gE,'dE,'^E] of the Euler system fll.ip . fll.2p . fl3.3p . fl3.4p if: 

• C is a bounded domain with a sufficiently smooth boundary, say dQ of class C°°-, 


• the initial data [po.Ejuo,e] satisfy 


go,E,'do,E £ lU^’^(f2), Uq^e £ i?^), Qq^e, 'do,E > 0 in fl; 


3 , 2 / 


(3.6) 


• the compatibility conditions 

afuo,E ■ n|an = 0 (3.7) 

hold for fc = 0,1, 2. 

We are ready to state the main result of this paper. 


7 












Theorem 3.1 Let VL <Z B? be a bounded domain with smooth boundary. Suppose that the thermo¬ 
dynamic functions p, e, and s are given by where pm, ^m, o^nd sm satisfy 

Let the transport coefficients p, 7] and X obey h2.10\) . h2.11\) . Let [qe,'&e,'^e\ be the classical solu¬ 
tion of the Euler system U.1H1.3\) . H3.3\) in a time interval {0,T), originating from the initial data 
uo,e] satisfying /\3.6\) . \3. 7p . Finally, let u] be a weak solution of the Navier-Stokes- 
Fourier system ([1.ITp, where the initial data satisfy 


go,'do > 0 a.a. in kl, J^qq dx> M, ||^o||L°°(n) + ||4IU°°(n) + ||uo||Loo(t^.^3) < D, (3.8) 


and where the scaling parameters a, v, u, and X are positive numbers. 

Then 

8 (^g,d,u gE,'dE,UE) (r) 

~ c(r, M, D) (^0, do, Uo go,E, do,E, uo,e) + max |a, z/, uj, X, ^ 

for a.a. r G (0, T). 


(3.9) 


a 


1 / 3 - 


V u^X, 


In view of the coercivity properties of the relative energy specihed in fl3.2p . Theorem 13 .1 1 provides 
an immediate corollary. 

Corollary 3.1 Under the hypotheses of Theorem \3. 1\ suppose that 

OJ 

a, z/, a;, A —)■ 0, and — 
a 


Then 


ess sup / 

re(o,r) 


V a 

—)■ 0, —)■ 0, __ ^ 0. 

ya yu^X 

(3.10) 

+ \q - + g\d - dEl] dx 

(3.11) 


< c{T, D, M)A (^a, z/,a;, A, ||po — ^o,E||L°°(r 2 ), H-do — 'do,£;||L°°(!rz), ||uo — , 

where A is an explicitly computable function of its arguments, 

A {a, z/, u), X, IIpo — ll-do — 'do,E||L°°(o), ||uo — uo,£;||loo(q./j3 )^ —)■ 0 


provided a,i',uj,X satisfy H3.10\) . and 






































Remark 3.1 As already mentioned above, the convergence stated in Corollarv \3.1\ is path dependent, 
the parameters a,h',u,\ are interrelated through It is easy to see that holds provided, 

for instance, 

12 3 

a —)■ 0, V = a", u = A = a'^, where /3>1, - <«<-, 0<7<1 — -a. 

^ O ^ 

Remark 3.2 The specific form of the function A is given in terms of the relative energy functional 
£ introduced at the beginning of this section. The constants in 113.9\} and \3.11\) depend also on the 
properties of the limit solution ug]. 

The rest of the paper is devoted to the proof of Theorem 13.1[ Possible extensions are discussed 
in Section [71 


4 Relative energy 

As shown in [10], any weak solution [g, il, u] to the Navier-Stokes-Fourier system with the boundary 
conditions fll.lOp . fll.lTp satishes the relative energy inequality 


£ (^,t9,u|r,0,u) 


t=T 

t=0 


(4.1) 


r r Q 

lo Jn d 


: V„u - 

< r 


V,l)' 


dx dt + A 


10 Jn 


Jo Jn 
Va;u) : Va;U dx dt — 


d j Jo Jn 

^(u - U) • V„U ■ (U - u) dx dt 


lur dx dt 


/o Jn 


d 


■ Va;0 dx dt + A 


u ■ U dx dt 


/o Jn 


10 Jn 


to Jn 


— • Va;0 dx dt 

+ U ■ ■ (U — u) dx dt — y j p{g, 'd)div 2 ;U dx dt 

d) — s(r, Q)^dtQ + g(^s{g, d) — s(r, 0))u • Va;0) dx dt 

+ Jo X ( (^ ~ f) f ^ ®)) 


for any trio of (smooth) test functions [r, 0, U] such that 

r, 0 > 0 in n, U • n|aQ = 0. 


(4.2) 
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Remark 4.1 Following Lions’ idea J77| / proposed for the incompressible Euler system, we may define 
dissipative solutions to the Navier-Stokes-Fourier system requiring solely certain regularity of [p, d, u] 
and the relative energy inequality to he satisfied for any trio [r, 0,U] as in Such an 

approach was used by Jessie, Jin, Novotny to attack problems on unbounded spatial domains. 


5 Energy estimates 


In what follows, we use the notation 


a ~ 6 if a < c6 


where c > 0 is a constant independent of the scaling parameters a, u, u, A, and of D and T. Similarly, 
we dehne a b and a ^ b. 

We start with some auxiliary estimates that follow directly from the structural hypotheses im¬ 
posed on the functions p, e, and s, see [9l Chapter 3] for the proofs: 


0 < esM(e,<?) ~ e(i + |iog(ri| + [iog()»)]+). 


(5.1) 

(5.2) 


Taking r > 0, 0 > 0 constant and U = 0 in the relative energy inequality fld.ip . and keeping 
flS.ip . fl5.2p in mind, we deduce the standard energy estimates: 


esssupig(o,T) In •) dx < c{D), 

esssuptg(o,T) In •) dx < c(D), 

esssupig(o,T) Iq Q^{t, •) dx < c(D), 

_ a esssup4g(o,(r) In •) dx < c(T>); _ 
together with the estimates following from dissipation and hypotheses fl2.10p . fl2.1ip : 


(5.3) 


^/o If! 


Va;U -I- V^u — |diV 3 ;Ul 


dx dt < c(D), A Jq /q |up dx dt < c{D), 


In + I log(^?)p] dx dt < c(B). 


(5.4) 


Now, by virtue of hypothesis (13.81) and the bounds (15.31 )i^ 2 , (I5.4p i . we may us a generalized 
version of Korn’s inequality P, Theorem 10.17] to obtain 


ll\/^u||j;^2(o,T;Wl.2(f7;_R3) < c{M, D)] 


(5.5) 
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whence, my means of the standard embedding relations for Sobolev functions, 


II V^u||i2(0,T;L6(O;H3)) < c{M, D) . 


(5.6) 


Moreover, in view of the interpolation inequality 

II II ^ II l|3/4 II ||l/4 


we obtain 

^,3/8;^1/8u 


< 


therefore, by virtue of (15.41 )i. 


Il\/j"u|| L6(Q;R3) IIVAuII L2(n;/J3) ~ (ll\/^^IU®(0;ij3) + ||\/Au||x,2(f^.^3)) , (5.7) 




,3/8 y 1/8 


U 


L2(o,r;L4(Q;ij3)) 


< c(M, D). 


(5.8) 


6 Convergence 

Our goal is to show Theorem 13.11 The obvious idea is to take 

r = qe, & = '&E, U = Us (6.1) 

as test functions in the relative energy inequality (14.ip . We emphasize again that such a step is 
conditioned by our choice of the slip boundary conditions for the velocity field in both the primitive 
and the target system. 

Next, we fix positive constants g, g, d, 'd in such a way that 

0<p< inf _Ps(t,x)< sup gE^tyX) <'g, 
te[o,r],a:er2 te[o,r],xGO 

0 <d< inf _'ds(t,x)< sup 'ds(t,x)<'d. 

te[o,T],xen te[o,r],xGr2 

Following [9] we introduce a decomposition of a measurable function F into its essential and residual 
part, specihcally. 


F = + m,,. IF]„ = HeE,'DE)F, |F],„ = (1 - <i(eE,'DE))F, 

where 

<h G (^^((O, cxo)^), 0<$<1, $ = lon the rectangle [p, p] x [i/, i?]. 
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Now, it follows from hypotheses 02.41) . 02.8p . and the estimates stated in 03.2p that 

lib - ^E]essllL2(o) + llb-^?E]esslli2(o) + II [u - ~ ^ ^9, u| Pfi, 7?^, U^; 

and, by virtue of 05.ip . 05.2p . 

I p|u —U£:Pdx+ / 1 + + p-d + a'd^ dx ~ fp, i?, u •de, u^;) . 

Jn JO L Jres V / 

Using the ansatz 06.ip we examine separately all integrals on the right-hand side of 
obviously, 

p(u - ue) ■ VxUe • (ue -u) dx dt < 8 (p, d, u qe, 'Oe, ue) , 

the relative energy inequality 04. ip can be written as 


t=T 

H=0 


^f§(^?,V.u):V.u 


In 1 ? 


8 (p,d,u Qe,i^e,Ue) 

qWV.tf).V^'| r /■ |„|2 <j,, 

/ JO Jn 




< 

r\j 


( 6 , 2 ) 

(6,3) 
. As, 


(6,4) 


10 Jn 


S('d, V3;u) ; V^Ue dx df — 


8 (^g,'d,u QE,'dE,UE) dt 

qi^, dx dt + A /" / u • Ue dx 

Jo Jn 


/o Jn iJ 


dt 


/o Jn 


^ lo In ~ ~ 

)(^dtUE + Ue ■ Vj;Ue) ■ (ue - u) dx df 


/o Jn 


p(g, i9)div^UE dx df 


0 Jn 


(p(s(p, iJ) - s(pe, ^e)) [dt-dE + Ue • Vx'Ue') dx dt 

f f ( (l - — ] dtp{gE,'dE) - —u-'\/xPigE,'dE)] dx dt. 
Jo Jn W gE J gE J 
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6.1 Integrals depending on viscosity and the heat flux 


On one hand, we have 




§(d, Va:U) : VxU ~ U 


2 ^ \ 
VxU + V^u --div^ul +-^|div^un, 


while, on the other hand, 
1/2 


= U 


§(d, V^u) : VxVLe 

(^xU + V^u - ^diVa;Ul^ : + u div,i;U div 


(6.5) 

( 6 . 6 ) 

xUe 


< du 


' m 

d 


Vj;U + V^u - -divj-ul 


+ 


rj^d) 

~r 


|div,j,un + z/c((5)d(^/i(i?) + 77 (d)) |Va 


Us 


< 


Id 


-^S(d, V^u) : V„u + i/c((5)d(/7(d) + r 7 (d)) | V, 


Us I 


provided the parameter h > 0 has been taken small enough. Moreover, in view of hypotheses fl2.10p . 
03.101) and the energy bound 05.31 h. 


V f +t;(!)))|Vj,UE/ di 


(6,7) 


< 

rsj 


Similarly, 


C>(s) + J v{P dx ~ 0 {y) + s^d"^ dx^ ~ C>(z/) + O 

dsq(d, v^d) ■ v^d > «(d) 2 > flvy 012 , ivy , /qM 2 \ laQ\ 

--~n;(^|V^d| +|V^log(d)| J, (6.8) 


while 


q(d,Vxd) ^qii1/2v7 q 

- ^ -= V^d-(fi:(d))^ Va;ds 

< ,,^fi:(d)|^^ 2 , /jcN /QMV 7 Q 12 / ldsq(d, V^d) • V^d 2 


(6.9) 


+ a;c(h)fi;(d)|V^ds| < -77- n 

2 V 


d 


+ n;c(h)fi)(d)|V^ds| 


where, furthermore. 


u [ h{'&)\'Vx'&e\‘^ dx ~ 0(uj) + — [ ad^ dx ~ 0(uj) + O f—^ 
JO a Jo Va/ 


( 6 . 10 ) 
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Finally, seeing that 


A / u ■ Ue dx < — [ |up dx + 0{X), 
Jn 2 Jn 


( 6 , 11 ) 


we may combine (16.5116.111) to reduce (16.dp to 
8 (^g,'d,u gE,'dE,UE'^ 


t=r 

t=0 


J 8 (^g, d, u gE^'dE, ue) dt + c{E, M )0 ^max |z/, a;, A, 

^ lo In ~ ~ 

[ [ g(dtUE + Ue ■ Vj;Ue) ■ (u^ - u) dx dt - [ [ p{g, ^9)div3.UE dx dt 
Jo Jn ^ ' Jo Jn 

-Jo J^{g{s{g,'d) - s{gE,dE))[dt'dE + ue -V^-dE) dx dt 

+ j f ((l- —) dtp{gE, ^e) - —u ■ V^p{gE, dx dt. 

Jo Jn \\ pe J Pe j 


( 6 . 12 ) 


6.2 Entropy integral 


The first integral on the right hand side of (16.121) is the most difficult one to handle. It can be written 
in the form ^ 

J^g[s{g,'d) - s{gE,'dE))(uE - uj ■Vx'dE dx dt (6.13) 

= Jo Jn^ “ ^^dE, ^?E)]ess {^E “ u) • V x'^E dx dt 

^ I In^ [■s(^, d) - s{gE, ^9i?)]res (ue - u) • Vx^E dx dt, 
where, in view of (16.21) . 

J^ J^g[s{g,'d) - s{gE,'dE)]ess{^E - ■'^x'dE dx dt 8 (^g,'d,u gE,'dE,UE) dt. 

Next, in accordance with (15.ip . 


J^ J^d [s(^>, id) - s{gE, ^9E)]res “ u) • dx dt 


(6.14) 


< 

r\j 


0 Jn 




(|u| + 1) dx dt + 


10 Jn 


1 + I log(^)| + [log(79)]+ v^|u - ueI dx dt, 
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where, by means of fl6.3p . 


/o Jn 


< 

r\j 


[l + I log(^)| + [log(t9)]+J^^^ v^|u - ue\ dx dt (6.15) 

/ / g + + Qi!^ dx df + / / ^|u — U£;P dx df ~ / S (g,'d,u gE,'&E,UE) dt. 

Jo Jn L Jres Jq Jq JO ^ ' 

Finally, by Holder’s inequality and fl6.3p . 

, 4/3 \ 3/4 

V 


r [ r \ ^ 


where, furthermore. 


,4/3 


dx = 


,1/3 


[ a-d'^ dx ~ c{E) 
Jn 


1/3 


Xo 1 jV 2AV6 iyl/2Al/6 

Summing up the previous estimates, we may rewrite fl6.12p as 


y/ 


0 . 


8 qe^'^e^^e) 


t=T 

t=0 


(6.16) 


< 

rs^ 


J 8 (^g, d, u gE, "de, dt + c{E, M)0 ^max 1 1 /, u, A, 


u UJ 

- 5 5 


1/3' 


/a a \y/i^J 

g{dtUE + ue • Vj,ue) • (ue - u) dx dt - ^ 

~ lo X ~dx dt 


g 


10 Jn 


1 - — ) dtp{gE,'dE) - —u • VxP{gE,'dE) 1 dx dt. 

. gE J gE J 


6.3 Remaining integrals 

Keeping in mind that [^EidEjUE] is a smooth solution of the Euler system, we may handle the 
remaining integrals on the right-hand side of fl6.16l) as follows: 


J J -l- ue • Va;UE) • (ue — u) dx dt — y J p(^,'d)div 2 .UE dx dt (6-17) 


lo Jn 


g{s{g, it) - s(^E, ^e)) dt^E + ue • V^di 


dx dt 
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Jo Jn 
Q 


1 —^ 1 dtp{gE, '&e) —■ VxP{gE, "&e) ] dx dt 

Qe J 


A. 

gE. 




—'^xPM{gE,'dE) ■ (ue -u) dx dt - [ f p{g,'d)diVxUE dx dt 
gE Jo Ju 

~ L X ~dx dt 

1 - — ) dtPMipE, ^e) - —U ■ VxPAiigE, 'Oe) ) dx dt + 0(a) 

gE) gE J 

= f f IpMigE^'dE) - p{g,'d)]diVxUE dx dt 

Jo Jn 

- [ [ \ l - —] [dtPM{gE,'dE) + ue ■VxPM{gE,'dE)] dx dt + 0{a) 

Jo Jn y gE J 



0 Jn 


- PM{g,'d)]^^^dwxUE dx dt 
■ [ [ g{sM{g, 'd) - SM{gE, ’^e)) [dt'^E + ue • Vx^?E]ess da: dt 



1 -^ 


JO Jn 


, [dtPAiigE, + Ue • '^xPuigE, ^i?)]ess ^a: dt 
gE / 

■djU ^Ej'^EjUe) dx + C>(a). 


write 


Jo Jn ^ ' 

Now, using the equations satisfied by [^Ej'^EiUe], together with Gibbs’ relation (ll.bjl . we may 
ite 

dtPM^gE, '^e) + Ue • VxPM{gE, '^e) 

dpAligEy'^E) /o , . OpMigE^'dE) /o Q 

--- [OtgE + Ue • VxgE) H-- [Ot'di 

og oxJ 


^E + Ue ■ 'Vx'&e) 


dpMigE,^E) ■,. 2 ^e) q , vy Q N 

gE --dlV,j,UE - gE -- [OtVE + Ue • Vx^e) ; 

og og 


whence (16.17^ gives rise to 


0 Jn 


[dtUE + Ue • Vo^ue) ■ (ue - u) dx dt - ^ 

-Jo J^g{s{g,'&) - s{gE,^E)) dt'&_ 


'd)diVa;UE dx dt (6.18) 


'e + Ue • Vx^E dx dt 
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< 

rsj 


< 

rsj 


+ dtp{gE,'>^E) - -^u-VxP{gE,'dE)^ dx dt 

f f PM^gE, ^e) - [qe - g) - PM{g, d) div^Ug dx dt 

Jo Jn op 

L J ess 

- / / g{sM{g, '&) - SM{gE, '^e)) [dt^E + dx dt 

+ ^ J^{g - [dt-dE + ue ■ dx dt 

(^g,'&,u gEy'dE^'ii-E) dx + 0{a) 



< 

r\^ 


div 3 ;UE dx dt 



- ess 


< 

rs^ 


Jo Jn 



./ U «/ i i 

f f PM{gE, ^e) - {qe - g) - PM{g, 

io Jn [ og 

gE{sM{g, 'd) - SM{gE, ^e)) [dt^E + Ue ■ Va;^?£]ess dx dt 
^ J^{g - gE)gE^^^^J^^^^ [dt^E + dx dt 

S (^g,d,u gE,'dE,UE'j dx + 0 {a) 

diVj;U£; dx dt 

ess 

hl^E + Ue ■ '^xdE]ess ^3^ dt 


0 Jn 


0 Jn 


f j puigE, ^e) - (gg - g) - puig, id) 

Jo Jn [ og 

r(sM(^,'d) - {g- _ SM{gE,d 

^ lo {d,d,u gE,'dE,UE) dx + 0{a) 

/ n N dpM{gE,'dE) , N dPuigE^^E) , ^ qs , qs 

PM{gE, dE) -- [gE — g) -- [ve — d) — PM{g, d) 

( { q\ ( xdsM{gE,dE) .q q ^dsM{gE,dE) 

gE[SM{g, d)-{g- gE) - ^ - [d- dE) - ^ - 


Jo Jn 



'o Jn 


dg dd 

X [dtdE + U^; • VxdE]ess ^3^ dt 

dpM{gE,dE) 


diVa-U^; dx dt 

ShligE, dE)'^ X 


dd 


-{dE-d) 


divj-u^; dx dt 


ess 
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0 Jn 


QE{'d - dx dt 


^ lo dx + 0 {a). 

Using the quadratic estimates fl6.2p . the above simplihes to 


/ n N dpM{QE,'dE) , X dpM{QE,'dE) , „ „x , qx 

Pm{Qe, Ve) -- [Qe — Q) - ^ - [ve — v) — Pm{q, '&) 


0 Jn 




dg dd 

dSuigE-, '&e) 


diVa;U£; dx dt 


(•d — df 


, dsM^QE, dfi) 


0 Jn 


dg ^ ad 

X [dt'^E + Ue ■ Va^d^Jess dt 
dpM{gE, dfi) 


SuiOE, 'dE)J X 


0 Jn 


PE(d - di? 


ad 

dsM^gE, d^ 

ad 


(dij - d) 


diVjxUE dx dt 


[dt^E + Uij • dx dt 


< 

rs-f 


^ J^S (^g,'&,u gE,'dE,viE) dx + 0 {a) 
dpniPE, d B) 


0 Jn 


ad 


-(ds-d) 


diVa;U£; dx dt 


0 Jn 


Finally, seeing that 


(di? - [dt{gEdE) + dw^{gE'dEVLE)] 


dsMigE^'^E) deM^gE'i'&E) / q x 

di?- ttt;-=-^ 77^-= C^(Pi7,dE), 


dx dt 


ad 


ad 


and, in accordance with fl3.5l) . 


Cv{gE,^E) [dtigE'dE) + div^{gE^EUE)] + ^e div^-UE = 0, 


ad 


we may infer that 


0 Jn 


](dtUE + Ue ■ VxU£;) ■ {ue - u) dx dt - 


, d)diVj;U£; dx dt 


(6.19) 
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da; dt 


~ lo ~ dt'&E + Ue -Vx^E 

+ f [ ((l -— ] dtp{gE,'&E) -—u-'\/xPigE,i^E)] dx dt 
Jo Jn \ \ gE J gE J 

S (^g^iJ^u gE,iJE,UE) dt + 0{a). 

Summing up the estimates fl6.16l) . fl6.19p and applying Gronwall’s lemma, we conclude 

S (p,t?,u (r) 


c{T, M, D) 8 (po, 4, uo go^E, 'do,E, uo,£;) + max a, u, oj, A, 


u 


a 


for a.a. r G (0,T). This relation, together with the bounds established 
flS.lip . We have proved Theorem I3.1l and Corollary I3.1I 


cu 


a 


1/3' 


a’ \Ve^ 

in 06.2P ' 


give rise to 


7 Concluding remarks 

The piece of information provided by 03.9p is slightly better and dehnitely more explicit in terms of 
the convergence rate than stated in Corollary 13.11 One could possibly do slightly better with respect 
to the values of the scaling parameters but it seems hard to remove the path dependence of the limit 
process. 

The same result with obvious modihcations in the proof can be obtained for a more general class 
of viscosity coefficients satisfying 

\p\'d)\ < c, p(l + < pdd) < 71(1 + 1 ?^), 0 < pdO) < rjil + 1 ?^) for all 'i9 > 0, with - < b < 1, 

— 5 

in place of fl2.10p . 

Similarly to the incompressible case, the situation described by the no-slip conditions 

u|ao = 0 

for the velocity in the Navier-Stokes-Fourier system remains an outstanding open problem. 
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